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Future Perspectives in Lattice Field Theory* 

C. T. Sachrajda, 

Dept. of Physics and Astronomy, University of Southampton, 
Southampton, S017 1BJ, UK 

I review some of the contributions which lattice simulations are likely to make during the next five years or 
so to the development of our understanding of particle physics. Particular emphasis is given to the evaluation 
of non-perturbative QCD effects in experimentally measurable amplitudes, and the corresponding extraction of 
fundamental parameters. 



1. Introduction 

I have been asked to talk about the future per- 
spectives of our field, which I find to be an ex- 
tremely daunting task. There are 388 partici- 
pants at this conference, each with his or her own 
view on how the subject may or should develop 
and all I can do is to present one of these 388 per- 
spectives. The focus of this talk will be on the role 
that lattice simulations arc playing now, and will 
continue to play in the future in the development 
of particle physics in general and in phenomenol- 
ogy in particular. The most important issue is 
how well it will be possible to quantify the non- 
perturbative strong interaction effects in experi- 
mentally measurable quantities. This is needed in 
order to be able to deduce fundamental theoreti- 
cal information from experimental measurements 
(e.g. in order to determine the CKM matrix ele- 
ments from experimental studies of weak decays) . 

Much of this talk will be based on the delibera- 
tions of the panel set up last December by ECFA 
(the European Committee for Future Accelera- 
tors) to advise it on the requirements for high- 
performance computing for lattice QCD in Eu- 
rope and I gratefully acknowledge my colleagues 
on the panel for their insights and help Q. The 
panel's terms of reference stated that "the main 
objective of this study is to assess the high per- 
formance computing resources which will be re- 
quired in the coming years by European physicists 
working in this field, and to review the scientific 
opportunities that these resources would open." 
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Most of the effort in improving the precision 
of lattice computations of physical quantities will 
be based on unquenching, the inclusion of light 
sea-quarks. A key ingredient in attempts to es- 
timate the precision achievable in future simula- 
tions is an understanding of the scaling behaviour 
of the algorithms used to generate full QCD con- 
figurations (as the quark masses and lattice size 
and spacing are varied). I start this talk with a 
brief review of some recent studies of this question 
(section ^). The following sections contain some 
of the main physics questions which will be stud- 
ied using lattice simulations during the coming 
years. The role of lattice computations in stan- 
dard model phenomenology will be considered in 
section [?[ and will include a discussion of a selec- 
tion of physical quantities which are already be- 
ing computed frequently and also of other quanti- 
ties for which lattice calculations are just begin- 
ning or for which we do not (yet?) understand 
how to perform the computations, even in princi- 
ple. Sections ||-|6| contain brief discussions of the 
prospects for lattice simulations of some quanti- 
ties in hadron spectroscopy, QCD thermodynam- 
ics and non-QCD physics. Finally, section ^ con- 
tains some closing comments. 

Perhaps the main excitement at this conference 
concerned developments in formulations of chiral 
fcrmions on the lattice, and we have had two very 
interesting plenary talks on this subject [Q. I will 
not discuss these developments because, although 
their impact is likely to be very significant, it is 
still too early to quantify the effect that the new 
formulations will have on the physics discussed 
below. 
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2. Computing Resources and Lattice Pa- 
rameters 

In order to gain some insights into the precision 
which will be achievable in the next five years or 
so, it is clearly necessary to try to forecast what 
computing resources might be available and this 
has been reviewed by Norman Christ 0] . Once it 
has been established, that we are aiming at com- 
puting power of terafiops, or tens or even hun- 
dreds of terafiops, it is also particularly impor- 
tant to have a good understanding of the scaling 
behaviour of the algorithms used to generate dy- 
namical quark configurations. This is necessary 
to determine the lattice sizes and spaces and the 
values of the quark masses with which we will be 
able to simulate. Our understanding of the scal- 
ing behaviour is currently far from complete, and 
a significant effort will have to be devoted to the 
study and development of algorithms. 

At this conference S. Giisken presented the re- 
sults of a study based on the performance of algo- 
rithms in the Monte-Carlo Runs (at a lattice spac- 
ing of 0.08 fm, with 2 flavours of Wilson fermions) 
of the SESAM/TxL collaboration §|. These au- 
thors extrapolate the observed behaviour of their 
algorithm to larger lattices and lighter quarks and 
an illustration of their conclusions is given in ta- 
ble |[ The table shows the estimated cost for gen- 
erating 100 independent configurations in Tflops- 
months. 

Some authors view the estimates in table as 
being too optimistic. The CP-PACS collabora- 
tion have recently tried to estimate the computer 
time required for a large-scale full QCD calcula- 
tion, with the quality of data comparable to that 
of the present quenched QCD study on the CP- 
PACS §]. They use a renormalisation group im- 
proved gauge action and two flavours of degen- 
erate quarks with the clover fermion action and 
estimate the required time to be of the order of 
100 TFlops-year. For example, they estimate that 
they would require 409 days on a 131 TFlops ma- 
chine (e.g. on a machine of 4096 PU's with a pro- 
cessing power of 32 GFlops/PU) in order to ob- 
tain 25,000 trajectories on a (3 fm) 3 spatial lattice 
with light quarks with masses down to 15 MeV, 
for which m^/m^ ~ 0.4 (realized by working on 
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0.40 
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66-141 


0.30 


8.2-15 




42 


152-332 



Table 1 



Estimated costs for generating 100 independent 
configurations in Tflops- months For a ^ 0.08 
the two bounds are given to reflect the uncer- 
tainty in the behaviour with the lattice spacing. 
The results from lattices which have an extent of 
5 or more pion correlation lengths are printed in 
bold type. 

a 48 3 x 96 lattice with a = 0.067 fm). 

A similar study is currently underway in the 
US, being made as part of the preparation for a 
proposal to the SSI (Scientific Simulation Initia- 
tive), and the preliminary conclusion is also that 
the Giisken et al. estimate is too optimistic by at 
least one order of magnitude Q. 

In table |2| I present three estimates of time 
required to generate 1000 independent configu- 
rations on a 3 fm lattice, with lattice spacing 
0.067 fm and m^/rrip — 0.4. This is essentially the 
CP-PACS scenario described above. Note that 
the estimate by Steve Sharpe is for two flavours 
of staggered fermions, and experience teaches us 
that the time required for Wilson-like fermions 
is at least a factor of 10 longer. This discussion 
demonstrates the importance of further investi- 
gations into the scaling behaviour of algorithms. 
The different groups should start by agreeing on 
the behaviour in the range of lattice parameters 
which they are able to simulate. There will nat- 
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Authors 



Estimate 



Action 



150 Tflops-Years RGI, Clover 
13 Tflops-Years Wilson 
2.5 Tflops-Years Staggered 



CP-PACS |j 

GLS 1 
Sharpe 
Table 2 

Three estimates of the time required to gener- 
ate about 1000 independent configurations on a 
3fm lattice, with lattice spacing 0.067 fm and 
m-x/mp = 0.4. 



urally remain uncertainties when extrapolations 
are made to larger lattices, smaller lattice spac- 
ings and lighter values for the quark masses. In 
spite of the uncertainties, and the obvious ex- 
pense in computing resources required to perform 
full QCD simulations with light quarks, the esti- 
mates above do indicate that the exciting goal of 
performing unquenched simulations at the same 
level of precision as current quenched ones within 
the next five years or so is a realistic one. 

Norman Christ has reviewed the status and 
prospects for machines for lattice simulations || . 
For the purposes of this talk I will simply assume 
that at least some of the major collaborations will 
have access to machines with O(lOTFlops) sus- 
tained power within the next five years or so. 

3. Lattice Phenomenology 

In this section I will discuss the prospects for 
quantifying non-perturbative QCD effects for a 
variety of physical quantities of phenomenologi- 
cal interest. For some of these quantities we un- 
derstand very well how to perform the compu- 
tations whereas for many others we do not. In 
the former category are many quantities which 
are well measured experimentally, and for which 
computations are used to demonstrate that lat- 
tice systematic uncertainties are under control. 
This then gives us confidence in the reliability 
of computations of unknown hadronic matrix ele- 
ments and other quantities needed to control non- 
perturbative QCD effects in physical processes. 

It has become conventional to present much 
standard model phenomenology in terms of the 
unitarity triangle. The unitarity of the CKM- 
matrix (written in the Wolfenstein parametrisa- 
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Figure 1. Unitarity Triangle corresponding to the 
relation in eq.(jj]). 
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leads to six unitarity relations, including the one 
most frequently considered which is shown in 
Fig. |l| (where p = p(l — X 2 /2) and fj = r/(l — 
A 2 /2)). Different physical processes give different 
loci for the position of the vertex A. In principle 
from the intersection of these curves we can deter- 
mine the position of A 7 of if more than two curves 
do not intersect at the same point, we would de- 
duce that effects of physics beyond the standard 
model are present. In practice, however, hadronic 
uncertainties are sufficiently large that we arrive 
at a region of allowed positions for A. Fig. || con- 
tains the allowed region deduced from one global 
analysis of current data ||. The joint aim of the 
theoretical and experimental communities is to 
reduce this region (see below). 

The remainder of this section is divided into 
two parts. In subsection 3.1 I discuss some 
quantities which are computed frequently in lat- 



tice simulations whilst in subsection 3.2 I con- 
sider others, for which lattice simulations are ei- 
ther just beginning or for which some conceptual 
progress is needed before such simulations will be 
possible. 
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Figure 2. Currently allowed values of the vertex 
A of the Unitarity Triangle || . The three regions 
correspond to 5%, 68% and 95% confidence levels. 

3.1. Frequently Computed Physical Quan- 
tities 

In table [| I present our preliminary estimates 
of the uncertainties in current and future results 
for a selection of physical quantities and matrix 
elements which are frequently computed in lattice 
simulations Q. The estimates in future calcula- 
tions are based on being able to generate 1000 
configurations at a lattice spacing of 0.08 fm with 
m^/nip = 0.4 and a sketch of how the improved 
precision will reduce the allowed region for A is 
presented in fig. |^ (without assuming any im- 
provement in the precision of experimental mea- 
surements). 

The mass of the strange quark, m s , is an inter- 
esting quantity with which to study the effects of 
quenching. It has been known for a considerable 
time that in quenched calculations one obtains 
a different value for m s depending on whether 
one uses the K- or $-mesons || (equivalently 
one cannot reproduce the physical value of the J- 
parameter, where J is basically the slope of rriy 
vs m 2 PS , where V {PS) represents vector (pseu- 
doscalar) mesons fLOA ). It will therefore be an im- 
portant check of our control of full QCD computa- 
tions to observe whether this problem is cured. At 
this conference the CP-PACS collaboration pre- 



quantity 


Present 


Future 




Uncertainty 


Uncertainty 


a s 


10% 


2% 




20-25% 


5% 


m b 


5% 


1% 


B K 


15% 


5% 


Other B's 


15-30% 


5% 


Id, Id s 


15% 


5% 


Jb, /s s , IbV b b 


20% 


7% 


F D ^ M (0) 


15-20% 


5% 


pB^M ( 0) 


20-35% 


7% 


Table 3 



Estimates of the present and future uncertain- 
ties in lattice calculations for a selection of phys- 
ical quantities, the strong coupling constant, 
quark masses, the B-parameter of K — K mix- 
ing and other B-parameters of four-quark op- 
erators, the decay constants of heavy mesons 
and heavy— dight form factors at zero momentum 
transfer. 

sented some preliminary results with 2 flavours of 
sea quarks, indicating that this may be the case 
(see the talk by R. Mawhinney Jll|), and it will 
be interesting to observe future developments. 

Among the recent interesting analyses in stan- 
dard model phenomenology, Parodi, Roudeau 
and Stocchi have used the measured values of 
|ejc | , |l«b/V c b|, Arrid and the bound on Am s to 
obtain p = 0.202 ± and fj = 0.340 ± .035 @ 
(see also the analysis in ref. constraining 
significantly the allowed positions of the vertex A 
of the unitarity triangle. These analyses rely on 
lattice computations of non-perturbative QCD ef- 
fects, such as those in Bk (for which the authors 
use B K = 0.86 ± 0.09) and /s d \f~B~B~ d (for which 
they take 210 t. 32 MeV), and of course the relia- 
bility of the analyses depends on that of the lat- 
tice computations. These authors, however, then 
perform an interesting exercise, by not including 
the lattice value of each of the hadronic parame- 
ters in turn and instead determining it from the 
analysis. In this way they find 

B K = 0.87 to.lt and ( 2 ) 
f Bd VB~B~ d - 223±13MeV. (3) 
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Figure 3. Scenario for the allowed region for the 
vertex A of the Unitarity Triangle using the esti- 
mates of future errors as in table |^ || . The three 
regions correspond to 5%, 68% and 95% confi- 
dence levels. 



errors) implies that direct lattice calculations only 
yield these form-factors at large values of the mo- 
mentum transfer (i.e. near the zero-recoil point). 
The extrapolation of the lattice results to smaller 
values of the momentum transfer leads to an im- 
portant source of systematic errors, even though 
we do have some theoretical guidance from the 
heavy quark effective theory (HQET), chiral per- 
turbation theory and from axiomatic properties of 
field theories (such as analyticity and unitarity) 
in performing these extrapolations [fig] . The er- 
rors due to the extrapolation will be significantly 
reduced by increased computing power (allowing 
us to go to smaller lattice spacings and hence to 
larger momenta for the mesons) . This is an illus- 
tration of the fact that for heavy flavour physics, 
increased computer power may also be fruitfully 
used to perform larger quenched computations. 
It is also to be hoped that there will be progress 
in understanding the scaling behaviour of form- 
factors at small momentum transfers to make the 
extrapolations more controllable. 



It is interesting to see the parameters which we 
are used to calculating in lattice simulations be- 
ing obtained in this way, albeit relying on lattice 
values for the remaining matrix elements and as- 
suming no new physics. 

One of the decay constants, Jd s1 has been mea- 
sured directly by several experiments, which the 
authors of ref. combine to obtain ^| 

f Ds = 241 ± 32 MeV , (4) 

in excellent agreement with lattice predictions, 
(e.g. in ref. jl| we quoted f Ds = 220±30 MeV as 
the average value from lattice simulations). The 
experimental result for frj s is therefore a signifi- 
cant encouragement to lattice phcnomcnologists. 

Finally in this subsection I consider the form 
factors for exclusive heavy — > light semileptonic 
and radiative decays, and in particular for B — > 
7r, p + leptons and B — > K*-f decays. The require- 
ment that the B and final-state light mesons have 
small momenta (in order to avoid discretisation 

2 Note that the Particle Data Group, more cautiously, 
present the spread of values from all experiments for 
in. @. 



3.2. Calculations which are Performed 
Less Frequently: 

The quantities considered in table [| represent 
an important set of phenomenologically impor- 
tant parameters, but there are also many other 
matrix elements and physical quantities for which 
we need to control the non-perturbative QCD ef- 
fects. I discuss a selection of these in this subsec- 
tion. 

The major difficulty in lattice phenomenology 
is that we have no general method for dealing 
with multihadronic states and final state interac- 
tions (see for example ref. p7[| ). We can com- 
pute amplitudes for processes with two particles 
at rest, which, when combined with chiral pertur- 
bation theory proves to be useful for kaon decays 
(but not for B-decays). I now consider exclusive 
nonlcptonic kaon and B-decays in turn. 

3.2.1. K — > 7T7T Decays: 

K — > 7T7T decays were considered in some de- 
tail during the early period of lattice phenomenol- 
ogy. As a result of the difficulties which were en- 
countered in the evaluation of the corresponding 
matrix elements and also because of the develop- 
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ment of heavy-quark physics the emphasis of lat- 
tice simulations changed towards quantities such 
as those in table ||. More recently the activity 
in kaon physics has increased again and at this 
conference was reviewed by Y. Kuramashi Jl9| . 
An important difficulty when studying kaon de- 
cays using Wilson-like quarks is to control the 
chiral behaviour (and the corresponding subtrac- 
tion of power divergences in many cases). The 
improvement in lattice technology and computing 
resources, together with theoretical developments 
such as non-perturbative renormalisation Jl8[ and 
new formulations for lattice fermions, imply that 
kaon decays should now become a major area 
of lattice phenomenology. The motivation for 
this is further underlined by the new measure- 
ments of e'/e from the KTEV and NA48 collabo- 
rations |2C[ , 

- = (28±3.0±2.8)1(T 4 [21] (5) 

- = (18±4.5±5.8)1(T 4 ||, (6) 

which are sufficiently large that one might hope 
to control the expected partial cancellations be- 
tween the matrix elements of the operators which 
are conventionally called Oq and 0$ []. 

A quantitative understanding of the AI =1/2 
rule in kaon decays is also very desirable and is 
an important benchmark for lattice QCD com- 
putations (note the interesting result in ref . p4| , 
which however has large errors bars, and, as 
is frequently the case for staggered fermions is 
plagued by very large perturbative corrections in 
the matching factors). 

3.2.2. B -> MM Decays: 

There is a flood of data becoming available 
on exclusive decays of B-mesons into two light- 
mesons, potentially giving fundamental informa- 
tion about the CKM Matrix and CP-violation. 
This flow of results will increase still further as 
the ^-factories and other approved experiments 

3 Since the conference the Riken-BNL-Columbia collabo- 
ration have presented results from a study using domain 
wall fermions in which they find surprisingly that these 
matrix elements have the same sign and obtain a negative 
result for e'/e (~ -(12.2 ± 6.8)10" 3 ) 



start taking data. Again the difficulty is in con- 
trolling the non-pertubative QCD effects. All the 
problems of non-leptonic kaon decays are present 
again, but in addition, the use of chiral pertur- 
bation theory is not applicable in £?-decays. New 
ideas are needed urgently and it may be that it 
is necessary to combine lattice calculations with 
models in order to make progress (see for exam- 
ple ref. p5|]). Any new ideas can also be tested 
on the huge amount of experimental data which 
already exists for nonleptonic charm decays. 

The difficulty in controlling hadronic effects 
in nonleptonic B-decays serves to underline the 
beauty of the one solid- gold process of B — > 
J/^K S , which is essentially free of these uncer- 
tainties and which will provide an accurate deter- 
mination of the angle (3 at the 6-factories []. 

3.2.3. B-Lifetimes: 

Lattice calculations of rates for inclusive non- 
leptonic decays are beginning, and are likely 
to make an important contribution to standard 
model phenomenology (see the review talk by 
Hashimoto To leading order in the heavy- 

quark mass all beauty hadrons are predicted to 
have the same lifetimes, and the corrections to 
this prediction can be calculated using an ex- 
pansion in inverse powers of the mass of the b- 
quark [S (see also refs J2j| and j3(J for reviews 
and references to the original literature) . For the 
ratio of lifetimes of the Af,-baryon and the B- 
meson the expansion leads to the prediction: 

li^l = l + 0-2% + O(±) , (7) 
T(B d ) m% 

where the on the right-hand side indicates that 
there are no operators of dimension 4 which can 
contribute and the -2% is an estimate of the 
contribution of dimension 5 operators obtained 
by comparing the spectroscopy of charmed and 
beauty mesons and baryons. In view of eq. (Q) 
the experimental measurement, r(A(,)/r(S c j) = 
0.78(4) is very puzzling. One possibility is that 
spectator effects (i.e. effects in which the light 
quark constituent of the hadrons participate di- 
rectly in the weak decay process), which in the 

4 The recent measurement from the CDF collaboration 
gives sin(2/3) = 0.79 §§. 
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heavy quark expansion appear at 0(1/ ml) but 
which have a phase-space enhancement |31J1 , are 
sufficiently large to account for the experimental 
result. This can be checked in lattice calcula- 
tions, by computing the matrix elements of the 
corresponding four-quark operators. A recent ex- 
ploratory calculation shows that spectator effects 
are indeed significant, and give a contribution to 
the right-hand side of eq.(@) of -(6-10)% US 
This calculation can readily be improved, and it 
is important to do so to learn about the practi- 
cability of using the heavy quark expansion for 
predictions of inclusive nonleptonic decays. 

A related quantity is Ar s /r s for which the first 
lattice results where presented at this conference 
by the Hiroshima group p4[ . 

3.2.4. Lightcone wavefunctions: 

Light-cone wave functions contain the non- 
perturbative QCD effects in hard exclusive pro- 
cesses (such as the form-factors at large mo- 
mentum transfers |55|,[56| or nonleptonic B- 
decays |37f ) and again lattice computations can 
make a significant contribution to the evaluation 
of these effects. To my mind, however, there has 
been insufficient effort in this field. In addition to 
calculating the moments of these wave-functions 
in the traditional way, it will also be interesting 
to explore and develop new suggestions for com- 
puting the wave-functions themselves (3^] . 

The main aim of this section was to stress the 
wide range of QCD phenomenology to which lat- 
tice computations can make a significant impact. 
Of course there are many examples in addition to 
those presented here, such as deep inelastic struc- 
ture functions |39| , proton decay amplitudes jl9| , 
shape functions for inclusive heavy— dight de- 
cays pofl , the electric dipole moment of the neu- 
tron etc. 



4. Hadronic Spectroscopy 

Hadronic spectroscopy continues to be a bench- 
mark of central importance for lattice computa- 
tions and the emphasis is now clearly on perform- 
ing unquenched calculations (see the comprehen- 
sive review by Bob Mahwinney |ll[]). In order to 
extrapolate the results obtained directly from lat- 



Quantity Uncertainty 2003-5 
m{q') - m(r]) 30MeV 
0++ Mixing 10% 
m(l~ + ) 20MeV 
m(H) 20MeV 

Table 4 

Guesstimates of the likely uncertainties which 
might be achieved for a number of key spectro- 
scopic quantities by 2003-5 or so. 



tice simulations to those corresponding to phys- 
ical u and ci-quark masses, guidance is needed 
from chiral perturbation theory. Estimates of 
how small the masses in simulations need to be 
in order for chiral perturbation theory to lead to 
reliable extrapolations are typically in the range 
corresponding to m^/mp ~ 0.3-0.4 (see ref. 41 
and references therein) and we have seen in sec- 
tion |^ that this is at the limit of what we might 
hope to achieve in the next few years. A particu- 
larly important milestone will have been reached 
when we are able to observe and control p — > 7T7t 
decays, which given the P-wave suppression of 
the decay-rates will take some time. 

In table ^ I have presented guesstimates of the 
relative precision which might be achieved for 
some key quantities in spectroscopy during the 
next three to five years or so Q . A significant ef- 
fort will be required to get a reasonable handle on 
all of these quantities. Samples of very large num- 
bers of configurations (10 4 or more) are required 
to get good signals for flavour singlet mesons. For 
the exotic 1 h state Doug Toussaint reminded us 
that the lattice result of 1.9 GeV should be com- 
pared with the experimental value of 1.4 GeV, so 
that we need to study mixing with meson-meson 
channels (42|. Note that the current theoreti- 
cal uncertainty on the mass of the 1 h meson 
is about 200 MeV. 

It will also be interesting and important to 
start to do "nuclear physics" with lattice simu- 
lations. These studies will presumably start with 
the simplest systems, such as the deuteron or the 
H-dibaryon (bound-state of two A-baryons). For 
these computations to be possible large volumes 
will be required, but, in addition, the fact that 
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the binding energies are small, 0(10 MeV), im- 
plies that again, high statistics will be needed. 

5. QCD Thermodynamics 

The commissioning of the RHIC accelerator at 
Brookhaven and the approval of the Alice experi- 
ment at CERN means that lattice thermodynam- 
ics will have a major phenomenological role to 
play as well as a theoretical one. At this con- 
ference Fritjof Karsch reminded us that the key 
features of the quenched theory are well under- 
stood (the existence of a first-order phase tran- 
sition, the equation of state, a critical tempera- 
ture of about 270 MeV) and that here also the at- 
tention is focused on simulations with dynamical 
quarks. The NUPECC Report on Computational 
Nuclear Physics concluded that in order to study 
QCD thermodynamics in simulations with light 
quarks, a lattice spacing of about 0.1 fm and a 
volume of 100 fm 3 , about 10 Tflops- years of com- 
puting effort are required. It appears that the 
effects of quenching are very significant, for ex- 
ample the critical temperature decreases substan- 
tially as quark loops are included (T c — > 170 MeV 
for two sea-quark flavours). In simulations with 
computing power in the 10 Teraflops range, it is 
expected that it will be possible to determine the 
critical temperature with a 5% accuracy. 

A schematic diagram of the expected phase 
structure of QCD in the (T, /i) plane is presented 
in figure |] (prepared by Simon Hands), where 
/i is the chemical potential. Verification that 
the structure is indeed as expected, and in par- 
ticular investigations of the fascinating scenario 
of a superconducting phase (or phases) of QCD 
at high /i (for a recent review and references to 
the original literature see ref. [Q), represents a 
major challenge for the lattice community (see 
the review by Shuryak at this conference [[44[). 
At non-zero /x, the action is complex, and con- 
ceptual progress is needed to develop reliable 
lattice calculations. It is, however, possible to 
study other quantum field theories with similar 
expected structures (but for which the action is 
real) and over the next few years we can look for- 
ward to the insights which these simulations will 
yield. 
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Figure 4. Schematic Phase Diagram for QCD. 



6. Non-QCD Physics 

In this talk I have concentrated on lattice QCD, 
but we should bear in mind the many application 
of lattice field theory to other areas of particle 
physics, some of which we have heard about at 
this conference. I will now briefly mention appli- 
cations to the electroweak phase transition and 
to supersymmetric gauge theories; other appli- 
cations include the study of phase structures in 
general (e.g. in QED), other applications in elec- 
troweak physics (strong Yukawa sectors, vacuum 
stability and mass bounds, heavy and strongly 
interacting Higgs bosons), and studies in quan- 
tum gravity. For many non-QCD applications 
the developments in new formulations of lattice 
fermions are likely to play a key role. 

6.1. The Electroweak Phase Transition 

The motivation for lattice studies of the elec- 
troweak phase transition is to gain an understand- 
ing of the baryon asymmetry of the universe and 
this subject has been comprehensively reviewed 
at this conference by Fodor 45 1 . For Higgs masses 
greater than about 60 GeV, infrared problems in- 
validate the use of perturbation theory, requir- 
ing the use of lattice simulations. The phase dia- 
gram determined from lattice simulations shows a 
line of lst-order phase transitions in the (T c , win) 
plane, terminating at a 2nd order end-point at 
mji — 72 GeV. In order for sphaleron produc- 
tion to be the mechanism for the generation of 
the observed baryon asymetry, as expected in the 
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standard model, there must be a strong first order 
transition, and hence the mass of the Higgs boson 
must be less than 72GeV, in contradiction with 
the experimental bound of 95 GeV from LEP ^. 
It would be reassuring to include the fermions ex- 
plicitly in these simulations (which would require 
major computing resources) but already these re- 
sults present significant evidence that physics be- 
yond the standard model is required. 

In the minimal supersymmetric standard 
model the bound is weaker, M H < 110 GeV, just 
above the experimental limit from LEP. 

6.2. Supersymmetry 

Much of the phenomenology of physics "be- 
yond the standard model" is based on supersym- 
metric field theories. These theories, with their 
high degree of symmetry, have extremely interest- 
ing non-perturbative phenomena (an outstanding 
example is Seiberg-Witten duality). The subse- 
quent developments have dramatically deepened 
our understanding of non-perturbative features of 
gauge theories, although it is not yet clear what 
the implications for non-supersymmetric theories 
such as QCD are. This makes the potential inves- 
tigations of supersymmetric theories using lattice 
simulations very exciting. 

Lattice studies of supersymmetric gauge theo- 
ries are a formidable challenge and it is natural 
to start with simpler field theories such as N = 1 
supersymmetric Yang-Mills theory and to investi- 
gate details of the expected non-perturbative fea- 
tures of confinement and spontaneous chiral sym- 
metry breaking. The computational requirements 
are similar to those for full QCD and current stud- 
ies are in their infancy (teaching us for exam- 
ple that the SU(2) theory has spontaneous dis- 
crete chiral symmetry breaking caused by a gaug- 
ino condensate). Simulations with small fermion 
masses are expensive (the supersymmetric limit 
has vanishing gluino mass), and computing power 
of the order of a Teraflop sustained would enable 
most of the important theoretical questions to be 
answered (mass spectrum, the SUSY potential, 

5 72 GeV is the position of the 2nd-order end-point. For a 
sufficiently strong transition the bound is even stronger, 
ran < 40 GeV, even more clearly in contradiction with 
experiment. 



Ward identities). The N = 2 SUSY theory has a 
considerably more complicated parameter space 
and so a large amount of effort will be required 
to explore its phase diagram and so no final re- 
sults can be expected soon. It remains a very 
important long-term goal. 

7. Conclusions 

It is my belief that in order for our commu- 
nity to continue to receive support in general, and 
for state-of-the-art computing facilities in partic- 
ular, we have to be perceived as having a defi- 
nite role in the development of particle physics. 
This is certainly the case at present and in this 
talk I have tried to illustrate the wide range of 
physical quantities and processes for which lat- 
tice results are being used extensively to quan- 
tify non-perturbative effects (see also the talk by 
Mangano |§0[). In many cases the leading obsta- 
cle to determing standard model parameters or 
other fundamental information from experimen- 
tal measurements, is not due to experimental dif- 
ficulties but to our inability to control the lattice 
systematic uncertainties as precisely as we would 
like. In discussing how our computations will im- 
prove in the future, it is easy to say that we will 
do the best that can be done with the resources 
which will be available. This is probably true, 
but in my view insufficient, and I believe that we 
also need to try and predict the "physics-reach" 
of the next generation (or two) of simulations. 
No doubt these predictions will have to be re- 
vised with time, but we need our theoretical and 
experimental colleagues to have a realistic pic- 
ture of our expectations, which are considerable 
but nevertheless limited. Together with my col- 
leagues on the ECFA panel we have tried to think 
about the prospects and strategies for the next 
five years, and some of our preliminary estimates 
can be found above. I hope that this will be a 
useful contribution to a debate within the com- 
munity about the future of lattice field theory. 
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